Infinite-range Ising ferromagnet in a time-dependent transverse magnetic field: 
quench and ac dynamics near the quantum critical point 
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We study an infinite range ferromagnetic Ising model in the presence of a transverse magnetic 
field which exhibits a quantum paramagnetic-ferromagnetic phase transition at a critical value of 
the transverse field. In the thermodynamic limit, the low-temperature properties of this model are 
dominated by the behavior of a single large classical spin governed by an anisotropic Hamiltonian. 
Using this property, we study the quench and AC dynamics of the model both numerically and 
analytically, and develop a correspondence between the classical phase space dynamics of a single 
spin and the quantum dynamics of the infinite-range ferromagnetic Ising model. In particular, we 
compare the behavior of the equal-time order parameter correlation function both near to and away 
from the quantum critical point in the presence of a quench or AC transverse field. We explicitly 
demonstrate that a clear signature of the quantum critical point can be obtained by studying the 
AC dynamics of the system even in the classical limit. We discuss possible realizations of our model 
in experimental systems. 

PACS numbers: 73.43.Nq, 05.70.Jk 



I. INTRODUCTION 

Quantum phase transitions have been studied exten- 
sively for several systems^. More often than not, simple 
prototype systems such as the Ising model in a transverse 
field or rotor models bring out many important charac- 
teristics of such transitions^. In most of these studies 
so far, only the equilibrium situation, where the system 
is taken adiabatically across the quantum critical point 
(QCP), has been addressed^. It is well known that in 
such cases, the thermodynamic properties of the system 
can be charted by determining the critical exponents and 
hence the universality class of the transitions* 2 ^. How- 
ever, no such general criteria exists for non-equilibrium 
dynamics of quantum critical systems. 

In recent years, there has been both theoretical and 
experimental progress in understanding the behavior of 
non-equilibrium dynamics near QCPs. On the exper- 
imental side, there have been several experiments on 
ultra-cold atoms in optical lattices which realize such 
non-equilibrium situations**. On the theoretical side, 
there have been studies of one-dimensional Ising mod- 
els in a transverse field and their variants^, driven quan- 
tum spin chains^, and correlation functions for quench 
dynamics close to a quantum critical poinli The main 
additional difficulty over the equilibrium case that one 
encounters in these theoretical studies arises from the fact 
that the description of non-equilibrium dynamics near a 
QCP necessitates the knowledge of all or at least the 
first few excited states in the Hilbert space of the many- 
body Hamiltonian. This precludes a detailed study of 
most systems which can be realized experimentally. The 
notable exceptions are situations where the systems can 
be described by either integrable models*'-, or can be 
analyzed numerically by exact diagonalization methods 



applicable to finite-size systems^. 

In this work, we study both the quench and the AC 
dynamics of an infinite-range ferromagnetic Ising model 
in the presence of a transverse field near its QCP. Such 
a model, as we show in this work, can be represented 
by a single large spin. This particular feature of the 
model allows us to analyze its non-equilibrium dynam- 
ics near the QCP. Further, the value of this single spin 
increases linearly with the system size so that the large 
spin and thermodynamic limits coincide for the model. 
This enables us to describe the system in the thermo- 
dynamic limit accurately using classical equation of mo- 
tions for a single spin. In particular, we can develop 
a correspondence between the classical phase space dy- 
namics of a single spin and the quantum dynamics of the 
infinite-range ferromagnetic Ising model in the presence 
of a transverse field. In view of the long range dipolc- 
dipole interactions involved in many systems with order- 
disorder transitions driven by tunneling fields, the study 
of quantum critical behavior of this model is not just 
a matter of theoretical curiosity, but can have applica- 
tion in ferroelectric^ like KH2PO4 and ferromagnets like 
Dyf^HsSO^O^O. Further this model can also can 
be realized in two-component Bose-Einstein condensates 
(BEC) of ultracold atomA 

The main results reported in this work are the follow- 
ing. First, we obtain the equilibrium phase diagram for 
the infinite range model and also obtain its collective ex- 
citations using a Holstein-Primakoff approach for both 
the ferromagnetic and the paramagnetic phase. Second, 
we study the quench dynamics of the equal-time order pa- 
rameter correlation function (EOC) and obtain analytical 
expressions for its long-time behaviors. Finally, we study 
the dynamics of the EOC in the presence of an weak ex- 
ternal AC transverse field; we show that the response of 
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the system in the presence of the AC field near the QCP 
involves multiple frequencies and hence appears noisy in 
the time domain. This behavior, which persists even in 
the classical large 5 limit of the model, is to be contrasted 
to that in both the paramagnetic and the ferromagnetic 
phases far from the critical point where the response of 
the system involves only a few frequencies and does not 
exhibit such noisy behavior. There are some studies on 
classical Ising models in the presence of oscillating lon- 
gitudinal fields- and some mean field and Monte Carlo 
studiesi2*ii in the presence of an oscillating transverse 
field. However, these studies are for completely different 
models and none of them investigates the true quantum 
dynamics. 

The organization of the paper is as follows. In Sec. 
PI we introduce the model, chart out its phase diagram 
and study its equilibrium properties. In Sec. IIIII where 
we discuss the quench dynamics of the system across the 
QCP both classically and quantum mechanically. In Sec. 
[j"Vl we study the AC dynamics of the model. This is 
followed in Sec. [V]by a discussion of possible realizations 
of our model in experimental systems and conclusions. 



II. EQUILIBRIUM PROPERTIES 

We consider a system of N spin-1/2 objects governed 
by the Hamiltonian 

- E s * s j - r E s " > w 



H 



N 



*<3 



where Sf = haf/2, a = x,y,z are respectively the x,y 
and z components of the spin-1/2 operator represented 
by the standard Pauli spin matrices a a . Here we assume 
that J > (ferromagnetic Ising interaction). This Hamil- 
tonian is invariant under the Zi symmetry Sf — ► Sf, 
S i ~ S i' and S i -> ~ S i ■ (The Z 2 symmetry would 
not be present if there was a longitudinal magnetic field 
coupling to Note that the model in Eq. JIJ dif- 

fers from the one studied in Ref. Il2l where the spins were 
taken to be living on two sub-lattices, with Ising inter- 
actions only between spins on different sub-lattices. We 
consider here the ferromagnetic case, for which the static 
quantum critical behavior is simple and easily derivable, 
while some significant features of its dynamic critical be- 
haviors are also analytically tractable. We take T > 
without loss of generality since we can always resort to 
the unitary transformation Sf — > —Sf, Sf — > —Sf and 
Sf — > Sf, which flips the sign of T but leaves J un- 
changed. Eq. can be written as 

H = - i ( Sf ot ) 2 
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and we have dropped a constant (J/2JV) J2i( S i) 2 = J / 8 
from the Hamiltonian in Eq. © . In the rest of this work, 
we shall use units h = 1. 



A. Mean Field Theory 

We begin with a mean field analysis of the thermo- 
dynamics of the model described by Eq. JTJ. Denoting 
the mean field value m = ^\ (Sf) /N, the Hamiltonian 
governing any one of the spins is given by 



JmSL 
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This is a two-state problem whose partition function can 
be found at any temperature T. If (3 = 1/ksT, we find 
that m must satisfy the self-consistent equation 
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This always has the trivial solution m = 0. In the limit 
of zero temperature, there is a non-trivial solution if 
T < J/2, with \m\ = (1/2)^/1 -4r 2 /J 2 ; the energy gap 
in that case is given by J/2. If T > J/2, we have m = 
and the gap is given by T — J/2. Hence there is a zero 
temperature phase transition at r c = J/2. The Zi sym- 
metry mentioned after Eq. is spontaneously broken 
and < Sf > becomes non-zero when one crosses from the 
paramagnetic phase at T > J/2 into the ferromagnetic 
phase T < J/2. 

In the plane of (fc^T/J, T/J), there is a ferromagnetic 
(FM) region in which the solution with m ^ has a lower 
free energy (the Z% symmetry is broken), and a param- 
agnetic (PM) region in which m = 0. The boundary 
between the two is obtained by taking the limit m —> 
in Eq. ©. This gives 2r/J = tanh(/3r/2), i.e., 



k B T 
J 



r 

j 



In 



_2T/J\ 



i - 2r/jy 



(6) 



The phase diagram is shown in Fig. ^ 



B. Holstein-Primakoff Approach 

We shall now find the low-energy spectrum of the 
Hamiltonian in Eq. (0. The form of (J2J shows that 
the total spin Sf ot = S(S + 1) is a good quantum num- 
ber; S can take any value from N/2 down to 1/2 or 0, 
depending on whether N is odd or even. For each value 
of S, the Hamiltonian describes a single spin; we shall 
see below that the low energy properties of the system 
are dominated by values of S lying close to N/2. The 
multiplicities D(S) for different values of S (not count- 
ing the trivial multiplicity of 25 + 1 due to different 
values of Sf ot ) are given by the following expressions: 
D(N/2) = N C = 1, D(N/2 - 1) = N d - N C =N-l, 
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general, if p is an integer satisfying 3 < p < N/2, then 
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FIG. 1: Phase diagram of the model in mean field theory. 
FM and PM denote ferromagnetic and paramagnetic regions 
respectively. 



For N — > oo and p « N/2, the leading term in D(N/2- 
p) is given by N p /p\ . 

We shall now find the ground state and low-lying ex- 
citations of Eq. (0 in a sector with a given value of the 
total spin S, assuming that T > J/2. To do this, we use 
the Holstein-Primakoff transformation, 
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a f ^2S - ata 



Hot ~r '"-'tot — " V ^ — , (8) 

where [a, a)] = 1. Assuming that S is close to N/2 which 
is a large number, we make the approximation of ignoring 
ala with respect to 2S in the last two equations in JHJl; 
this gives 
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Slot = yf (a + T). (9) 
The Hamiltonian in Eq. J2J then takes the form 

if = - TS + T at a - 4^ (a* + a) 2 . (10) 
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The spectrum of this Hamiltonian can be found by trans- 
forming to the variables q and p, where 



ip 



and 
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vp 



V2 ' V2 
and [q,p] = i. Then Eq. (|1C)|> takes the form 
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This describes a simple harmonic oscillator with the fre- 
quency 
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[Since the maximum value of S is N/2, we see that us 
is real since T > J/2. If T < J/2, we have to analyze 
the problem differently as discussed below.] The energy 
spectrum of Eq. I|12f> is given by 
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where n = 0, 1, 2, • • • . The ground state of the system 
corresponds to S = N/2 and n = 0. 

Since 5 takes the values N/2, N/2-l, N/2-2, • • • , Eq. 
(|14f> shows that there are two towers of equally spaced 
excitations: one with a spacing of T coming from the 
first term (arising from excitations in which the total 
spin changes) and the other with a spacing coming 
from the second term (arising from excitations in which 
the total spin does not change) . Only the energy spacing 
r turns out to be thermodynamically significant; this is 
due to the multiplicities given in Eq. (Q as we will now 
see. The partition function of the oscillator in Eq. i|14|) 
at an inverse temperature /3 is given by 
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The complete partition function is therefore given by 



(15) 
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In the thermodynamic limit N — > oo, the sum in Eq. 
(|16fl will be dominated by values of S lying close to N/2. 
We can therefore write the partition function as 
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where we have approximated ui by its value at S — N/2, 
namely, co ~ r-\/l — J/2T. The free energy per spin is 
then given by 
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to lowest order in the quantity e~~P T , in the limit of zero 
temperature. This shows that the thermodynamic gap is 
given by T. Note that this agrees with the gap obtained in 
mean field theory. The harmonic oscillator energy levels 
with spacing w§ correspond to a collective excitation of 
all the spins and therefore they do not appear in mean 
field theory (which only takes into account excitations at 
a single site). 

We will end with a brief discussion of the case T < 
J/2. A classical analysis with S = {S x , S y , Sf ot ) = 
S (cos c/)sin 9, sin (f>sm9, cos 9), where 9 and (j> are the 
usual polar and azimuthal angles, gives the Hamiltonian 
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cos 2 9 — TS sin 9 cos < 
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The lowest energy configuration of this is given by S = 
N/2, 9 = wa~ 1 {2T/J) and = 0. There are two solu- 
tions for 9, one lying in the range [0, tt/2] and the other 
in the range [tt/2 ) tt]. We choose either one of these as 
the ground state and perform a Holstein-Primakoff trans- 
formation around it. We then find that there are again 
two towers of equally spaced excitations, one with spac- 
ing J/2 and the other with spacing (J/2)^/l - AT 2 / J 2 
corresponding to a collective excitation. The thermody- 
namically significant spacing is J/2. 

Note that as T approaches J/2 from either above or 
below, the collective excitation energy softens and goes 
to zero. Thus it is the collective excitation gap (not the 
thermodynamically significant gap) which is sensitive to 
the QCP lying at T c = J/2. 

We would like to emphasize here that our model has 
infinite range interactions and we are interested in the 
thermodynamic limit N — » oo. In that limit, the effec- 
tive Hamiltonian seen by any one spin is just the mean 
field produced by all the other spins. Hence there is no 
difference between the phase diagrams obtained by the 
mean field theory and the Holstein-Primakoff approach. 

A difference between mean field theory and the 
Holstein-Primakoff approach is that only the latter cor- 
rectly describes the collective excitation which involves 
all the spins (as we have discussed above). However, a 
single mode like the collective excitation does not play 
any role in determining the phase diagram in the ther- 
modynamic limit. 



III. QUENCH DYNAMICS 

In this section, we will study the quench dynamics of 
the spin model across the QCP. We saw above that it is 
the collective excitations (in which the total spin does not 
change) which are sensitive to the presence of a QCP. We 
will therefore restrict our attention to the lowest energy 
sector in which the total spin is given by S = N/2. 

A quench dynamics across a QCP has been studied ear- 
lier for Ising-like systems with short-range interactions 
and driven XY spin chains^. Here, taking advantage of 
the fact that the large S (and hence classical) limit of 
our model is also the thermodynamic limit, we shall de- 
rive an analytical expression for the long-time average 
of the EOC in the thermodynamic limit as the system 
is quenched across the QCP starting from the paramag- 
netic phase. In this section, we shall mostly follow the 
method of Ref. @. 



A. Quantum Analysis 

To begin with, we study the dynamics of the EOC (de- 
fined as ((St ot ) 2 ) /S 2 ) by changing the transverse field V 
from an initial value Ti > r c to a final value T/ sud- 
denly, so that the ground state of the system has no time 
to change during the quench. In this case, just after the 



quench, the ground state of the system can be expressed, 
in terms of the eigenstates \n) of the new Hamiltonian 
n f = -(J/4S)(Si ot f-T f Sf ot as 

|V) = ]>>», (20) 

n 

where c n denotes the overlap of the eigenstate \n) with 
the old ground state As the state of the system 

evolves, it is given at time t by 

\m) = 5>e-«^/A| n)) (21) 

n 

where E n = (n\Hf \n) are the energy eigenvalues of the 
Hamiltonian Hf. The EOC can thus be written as 

m)\ (Stotf/s 2 \m) = E c » c ™ cos - E ™) 

m,,n 

x(m\(S* ot ) 2 /S 2 \n).(22) 

Eq. (1221 can be solved numerically to obtain the time 
evolution of the EOC. However, certain qualitative fea- 
tures of the dynamics can be extracted from Eq. 1221) 
without resorting to numerics. First, we note that 
when Tf and Ti both lie in the paramagnetic phase, 
we expect the old ground state to have a large over- 
lap with the new one, so that c n ~ S n i. In this 
case, ({Sf ot ) 2 (t)) is expected to undergo small oscil- 
lations about ((Sf 0t ) 2 (t = 0)). On the other hand, if 
Tf <C r c <C T.^ there is very little overlap between 
the two ground states and the amplitude of oscillation 
is again expected to be small. Note that the situation is 
completely different from the adiabatic turning on of T, 
where for Tf <C r c , the system has a maximal value of 
((Stot) 2 )- I n between these two regimes, when Tf ~T C , 
the old ground state is expected to have significant over- 
lap with many eigenstates |m), and the oscillation ampli- 
tude can be expected to be large. These qualitative ex- 
pectations are verified in Fig. [21 Here, we have quenched 
the transverse fields to Tf/J = 0.9,0.01, and 0.4 start- 
ing from Ti/ J = 2.0. The oscillation amplitudes of the 
EOC for S = 100, as shown in Fig. |2 are small for 
Tf = 0.9 and 0.01, whereas it is substantially larger for 
T f = 0.4. 

Next, to understand the dynamics of the EOC in a 
little more detail, we study its long-time averaged value 
given by 

0= YimJ({Sl t ) 2 {t))) T /S 2 

= ^E c 'H(S t 2 ot ) 2 H (23) 

n 

for different T f . Note that the long-time average depends 
on the product of the overlap of the state \n) with the 
old ground state and the expectation of {Sf 0t ) 2 in that 
state. From our earlier discussion regarding the dynamics 
of ((Sjot) 2 ), we therefore expect O to have a peak some- 
where near the critical point on the ordered side where 
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FIG. 2: Dynamics of ({S? ot ) 2 ) / S 2 for S = 100 after quench- 
ing the transverse field to different values Tf/J from an initial 
field Ti/J = 2. The oscillation amplitudes are small, as seen 
from the solid (red) and dotted (blue) curves corresponding 
to Tf/J = 0.9 and 0.01 respectively, far away from the criti- 
cal point T c /J = 0.5. The oscillation is large in the ordered 
phase near the critical point as seen from the dashed (black) 
curve T s /J = 0.4. 



FIG. 3: Plot of the long-time average O as a function of Tf / J 
for different S. In the plot the solid (blue), dotted (black), 
dash-dotted (green) and the dashed (red) lines represents re- 
spectively the results for S = 50, 5" = 100, 5" = 200 and 
S — 500. O peaks around Tf/J = 0.25 and the peak value 
decreases with increasing S. For all plots we have chosen 
Til J = 2. 



such an overlap is maximized. This is verified by explicit 
numerical computation of Eq. (|23[1 in Fig. [21 for several 
values of S and Ti/J = 2. We find that O peaks around 
Tf/J — 0.25, and the peak height decreases slowly with 
increasing S. 



B. Classical Analysis 

To understand the position and the system size depen- 
dence of the peak in O, we now look at the thermody- 
namic (large system size) limit; in the present model, this 
is also the large S and therefore classical limit. With this 
observation, we study the classical equations of motion 
for S = S (cos (j>sm9, sin</>sin 9, cos 9) for T = Tf. In 
the present model, S is a constant. Thus in the classical 
limit, we need to study the equations of motion for 9 and 
<f>. To this end, we note that the classical Lagrangian can 
be written in terms of 9 and 6 &si& 



L 



-S fl - cosf 



dt 



H[6,<f>] , (24) 



This gives the equations of motions 
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cos6» + Tf cot 8 cos <j> . (25) 



Eq. I|25|) has to be supplemented with the initial con- 
dition that Sf ot =Satt = Q. The condition Sf ot = S 



corresponds to# = 7r/2, (f> = which is also a fixed 
point of Eq. 1)25(1. Therefore we shall start from an 
initial condition which is very close to the fixed point: 
9 = tt/2 — e, 4> = e, where e is an arbitrarily small con- 
stant. Further, since the motion occurs on a constant 
energy surface after the quench has taken place, we have 
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smf cos < 



(26) 



Using Eqs. 1251 and (|26|l . we get an equation of motion 
for 9 in closed form, 



d9 
lit 



T 2 f sin 2 (0) - [r 



j cos 2 9 



f(0). (27) 



It can be seen that the motion of 9 is oscillatory and has 
classical turning points at Q\ = sin -1 (|1 — 4Ty/J|) and 



c 2 = tt/2. One can now obtain {{Sf ot ) 2 ) T 
fromEq. $Tf$, 
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cos" 



<cos 2 0) = M/V 



V = 
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cos 2 9 



f(0) 

yj8T f (J - 2T f )/J 
1 



,19 



f(0) 



(28) 

(29) 
(30) 



When trying to evaluate T>, we find that the integral has 
an end-point singularity at 9i ; this can be regulated by a 



6 



cut-off rj so that 62 
V - 



7r/2 — rj. With this regularization, 
The cut-off used here 



-Jln(ry)/ v /r / (J-2r / )/2 
has a physical meaning and is not arbitrary. To see this, 
note that the angles (6, </>) define the surface of a unit 
sphere of area An. This surface, for a system with spin 
5, is also the phase space which has 25+1 quantum 
mechanical states. For large 5, the area of the surface 
occupied by each quantum mechanical state is therefore 
47r/(25 + 1) ~ 2tt/S. In other words, each quantum 
mechanical state will have a linear dimension of order 
l/\/~S\ this is how close we can get to a given point on 
the surface of the sphere. Note that this closeness is 
determined purely by quantum fluctuation and vanishes 
for 5 — > 00. Thus 77, which is also a measure of how close 
to the point = n/2 we can get, must be of the order 
of 1/ %/5 ; this determines the system-size dependence of 
(cos 2 6>) T . Using Eq. (J2HJ, we finally get 



(cos 2 9) rj 



i6iy (J-21V) 

J 2 ln(5) 



(31) 



Eq. (|31[) is one of the main results of this work. It demon- 
strates that the long-time average of the EOC must be 
peaked at Tf/J = 0.25 which agrees perfectly with the 
exact quantum mechanical numerical analysis leading to 
Fig. Moreover, it provides an analytical understand- 
ing of the 5 (and hence system-size) dependence of the 
peak values of Tf/J. A plot of the peak height of O as a 
function of l/ln(5) indeed fits a straight line, as shown 
in Fig. 21 So we conclude that the peak in O vanishes 
logarithmically with system size 5. 
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FIG. 5: Plot of the DFT power spectrum of the EOC 
{ip{t)\{St ot ) 2 \ip{t)) / S 2 for different values of V d . Panels (a) 
and (c) show the DFTs in the ferromagnetic and paramag- 
netic phases respectively, whereas panel (b) shows the same 
for the quantum critical region. In all three panels, the height 
of the peak at the drive frequency (u) = 0.1) is maximum and 
is scaled to 1 (chopped off in the figure). The plots clearly 
indicate contributions from a relatively large number of fre- 
quencies to the DFT power spectrum near the QCP. 



Sec. IIIII we shall first study the quantum dynamics of 
the system in the presence of a transverse ac field. This 
will be followed by an analysis of the classical equations 
of motion (Eq. Q25fl) which describe the behavior of the 
system in the large 5 limit. 



FIG. 4: Plot of the maximum peak height O max of the long- 
time average of the EOC as a function of l/ln(5). The 
straight line shows the linear fit. 



IV. AC DYNAMICS 

In this section we study the dynamics of the system 
in the presence of an oscillatory transverse field. As in 



A. Quantum Analysis 

In the presence of a transverse AC field, the total 
Hamiltonian is 



H{t) 



J 

~2N 



(S? ot ) 2 -[T d + T a co S (ut)]S? ot , (32) 



where denotes the static part of the transverse field, 
r a is the amplitude of the AC part, and u> is the drive 
frequency. Note that for small r a , Td determines the 
point in the phase diagram about which the AC dynam- 
ics takes place. For the rest of this work we shall restrict 
ourselves to this limit and for definiteness take a small 
ratio T a /Td = 0.18. We have checked that the qualita- 
tive picture of the AC dynamics that we describe below 
remains the same for a large range of this ratio as long 

as r a /r d < 1. 

To elucidate the properties of the AC dynamics of our 
model, we calculate the EOC (i(j(t)\(Sf ot ) 2 \ip(t))/S 2 . For 
a sufficiently small driving frequency u>, we are in the adi- 
abatic regime and the system always stays in the instan- 
taneous ground state of H{t). In this case, the behavior 
of the EOC is the same as that for the static case, and 
has been well studied^. However, the situation changes 
when we drive the system with moderately high frequen- 
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cies due to possible Landau-Zener type transitions be- 
tween different time-dependent levels of H(t), and a full 
out-of-equilibrium analysis of the problem becomes nec- 
essary. If, however, the driving frequency becomes very 
high compared to other natural frequency scales of the 
system, the eigenstates of H{t) vary much faster com- 
pared to the characteristic time scale of most of the pos- 
sible Landau-Zener transitions. Hence such transitions 
fail to take place and the response is characterized by 
a few trivial frequencies. In this work, we shall there- 
fore investigate the case of moderately high frequency 
and study the characteristics of the frequency spectrum 
of the EOC. For the sake of defmiteness, we shall choose 
uj = 0.1J in the rest of this work. 

For a given set of parameters {Td,T a and u>), we start 
with the ground state oiH{t = 0) (in Eq. as the ini- 

tial state \ip(t — 0)). We then solve the time-dependent 
Schrodinger equation to obtain the time evolution of the 
state \ip(t)), using 4th order Runge-Kutta method with 
adaptive step-size control^. The time variation of the 
EOC is found, and the discrete Fourier transform (DFT) 
of the data is taken to obtain its frequency spectrum. 
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FIG. 6: The time evolution (in units of ft/ J) of the overlaps of 
the instantaneous state \ip(t)) with the time-dependent eigen- 
states \n(t)) of the total Hamiltonian Hit) (n = 1 denotes the 
ground state, n — 2 the first excited state and so on). Results 
are shown for n = 1, 3, 5 and 7. 



The DFT power spectra show that inside the para- 
magnetic and ferromagnetic phases, the response of the 
system is dominated by a few frequencies (panels (a) and 
(c) of Fig. [SJ. These are either the driving frequency 
lu = 0.1J or the natural oscillation frequencies of the 
system as will be shown in more details in Sec. HVBl 
In contrast, near the QCP (panel (b) of Fig. ISJ, several 
frequencies appear in the spectrum. Thus we find that 
the EOC displays a noisy behavior in the time domain 
near the QPT. 

To look into the origin of these frequencies, we have 
plotted the instantaneous overlaps of the time dependent 
state \ip(t)) with the first few instantaneous eigenstates 



\n(t))'s of H(t), for to = 0.1 in Fig. These overlaps are 
significant only for n = 1 (ground state), n — 3, 5 and 7 
for all times. The overlap is negligible for all other odd 
n and vanishes for all even n by symmetry. (The states 
with odd and even values of n have opposite values of 
the Z2 parity mentioned after Eq. (JIJ. Hence the ma- 
trix element of S^ ot between states with odd and even n 
vanishes). In addition, we have also plotted the instan- 
taneous gaps G mn {t) = E n (t) — E m (t) in Fig. In the 
ferromagnetic region, as shown in panel (a) of Fig. 
the transitions are mainly between the two levels, n — 1 
and 3. However, the overlaps have much higher average 
value for n = 1 than that for n = 3. This indicates that 
the system mostly stays close to the ground state. Also, 
as can be seen from panel (a) of Fig. [7] the relevant 
gaps G mn (t) are almost constant, resulting in a periodic 
behavior dominated by the driving frequency u) = 0.1J 
and a frequency 0.52 J which is close to the almost static 
value of the gaps. Deep inside the paramagnetic region 
(panel (c) of Fig. EJ, the system is found to be in the 
ground state (\n — 1(£))) throughout the evolution; all 
other overlaps tend to zero. This is because in this phase, 
as can be seen from panel (c) of Fig0 the instantaneous 
gaps G mn (t) = E n (t) — E m (t) always remain large and 
prohibits occurrence of any transition from the ground 
state to the higher levels. Thus the dynamics is again 
governed by a few frequencies and has a simple periodic 
behavior. In contrast, near the QCP (panel (b) of Fig. |5J) 
all of the levels n — 1,3 and 5 have substantial overlaps 
with the instantaneous ground state (n — 7 also has a 
small contribution). The overlaps are much flatter and 
more evenly distributed among the levels n = 1 and 3 on 
an average. Further, as can be seen from panel (b) of Fig. 
[7| the gaps Gmn undergo appreciable oscillations, creat- 
ing possibilities of Landau-Zener transitions at multiple 
frequencies. Thus, quite generally, the AC dynamics near 
the QCP involves multiple frequencies and is expected to 
exhibit noisy behavior in the time domain. 

Before ending this section, we would like to point out 
that the time evolution of ((S*^) 2 ), for some specific 
parameter values, may have large oscillations with an 
unusually long time period and discuss the reason for 
its occurrence. To this end, we present an example of 
this phenomenon in Fig. |SJ The figure shows oscilla- 
tions with a time period of about 1100 which is much 
longer than the period of the oscillating transverse field, 
namely, 2?:/ to ~ 63. This phenomenon is well-known in 
the context of Rabi oscillations which can occur when a 
two-state system is subjected to a periodic potential 15 . 
For instance, consider a Hamiltonian which is a 2 x 2 ma- 
trix of the form (u>o/2)a z + (b cos u>t)a x . If b is small and 
we start initially in the ground state of the term (u>q/2)o z 
(the energy gap corresponding to this term is given by 
u>o), we find that if u is close to uj , the system period- 
ically makes transitions to the excited state and returns 
to the ground state. The time period of this oscillation is 
given by 2%/ — w ) 2 + b 2 ; this is much longer than 
2-7T /lu if lo/loq is close to 1 and b <C lu. The amplitude 
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FIG. 7: Plot of the relevant gaps G mn (t) = E n (t) - E m (t) 
between different eigenvalues of H(t) as a function of time 
(in units of h/J) for different values of IV In (a), in the 
ferromagnetic phase (T^ = 0.1), the gaps are almost static and 
equal at all times. In (b), near criticality (Fd = 0.55), the gaps 
vary considerably with time and attain quite low stationary 
values. In (c), in the paramagnetic phase (Yd = 1-0), the gaps 
vary considerably with time, but are always much larger than 
their counterparts in the other two phases. 



of the oscillations is maximum for u> = ujq. (Similar 
long time period oscillations can occur if lo/ujq is close 
to 1/2, 1/3, • • • ). Returning to our problem, for the pa- 
rameters given in Fig. [SJ we find that the system is quite 
well described if we only keep the first three eigenstates 
corresponding to n — 1, 3, 5 of the time-independent part 
of the Hamiltonian in Eq. (|32[l . If we truncate the Hamil- 
tonian to only these three states and study the AC dy- 
namics starting initially with the first eigenstate, we find 
that the probabilities of being in the three states have 
large and long time period oscillations as shown in Fig. 
El We thus have a three-state version of Rabi oscilla- 
tions. [One of the parameters used in Fig. [SJ r a , has 
been chosen to be different compared to the earlier calcu- 
lations. This is because the presence of Rabi oscillations 
is sensitively dependent on the parameter values; it was 
necessary to change the value of r a , keeping the other 
parameters fixed, to show a striking example of such os- 
cillations.] Note that although the probability of being in 
the state n = 5 is small at all times as shown in Fig. [5] it 
is necessary to keep this state in order to obtain the long- 
time oscillations. We have checked numerically that if we 
truncate the Hamiltonian to only the states n = 1 and 
3, we do not get these large and long-time oscillations. 
We also note that the driving frequency u> = 0.1 is close 
to one-fourth of the gap between the states n = 1 and 
3, namely, G13 = 0.408. The Rabi oscillations are found 
to peak at lu ~ 0.101 and to disappear if u is changed 
from that value by just 4%. Therefore it is expected that 
the low-frequency peaks in the DFT spectrum are not 
generic features at all drive frequencies, but will be seen 
only for special values of the drive frequencies. 
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FIG. 8: Plot of the DFT power spectrum in panel (a) and the 
time evolution in panel (b) of the EOC {4>(t)\(S? ot f\ij(t)} / S 2 
for S = 50, r d = 0.55, T a = 0.055, and u = 0.1. Note 
the oscillations with a very large time period of about 1100 
in panel (b), and the corresponding peak at a frequency of 
about 0.006 in panel (a). 
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FIG. 9: Time evolution of the probabilities of being in differ- 
ent states when the model considered in Fig. |H|is truncated 
to the three states n = 1, 3 and 5. 



B. Classical Analysis 

The classical analysis of the AC dynamics of the model 
involves a numerical solution of Eq. i|25|) with — > 
I'd + r a cos(wi). As in the case of the quantum analysis, 
we are going to restrict our analysis to an intermediate 
drive frequency lo = 0.1J and S> r a . Note that we are 
interested in the classical limit of the quantum dynamics, 
for which the system is in it's ground state at t = Q. Thus 
we choose the initial condition for the classical dynamics 
such that the system is close to its fixed point. These 
fixed points, as can be easily seen from Eq. (|25|l . are 
given by (9f,<fif) — (ir/2,0) when the system is in the 
paramagnetic phase with Yd > 0.5, and by {8f,4>f) — 
(sin _1 (2Fd), 0) for the ferromagnetic phase where Td < 
0.5. Note that the fixed point structure of Eq. 125|) 
changes at Td/J = 0.5 which is also the location of the 
QCP. Therefore, we choose the initial condition (8f — e, e) 
for small e = 0.1. We have checked that smaller values 
of e do not lead to qualitative change in the essential 
features of the dynamics. Also, in the rest of this section, 
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we shall maintain the same ratio T a /Td = 0.18 as in the 
previous section. 

Before carrying out the numerical simulations, it is use- 
ful to obtain a qualitative understanding of the dynamics 
of the system in the presence of the driving field. Since 
the amplitude of the drive field T a is small and the mo- 
tion starts from near the classical fixed point, we expect 
that a small amplitude analysis for the motion of 9 and 4> 
would give us a qualitative understanding of the dynam- 
ics for both the paramagnetic and ferromagnetic phases. 
Of course, this expectation is not met near the quantum 
critical point, and we are going to discuss this point in 
more detail subsequently. 
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and from the drive frequency u> — 0.1J. This expecta- 
tion can be verified from panel (c) of Fig. The exact 
location of all the frequencies, however, cannot be eas- 
ily deduced from this qualitative argument. Next, let us 
consider the ferromagnetic phase at T c i = 0.1. Here the 
fixed point is given by (0f,(f)f) = (sin - (2P//J),0). A 
similar linearization about the fixed point leads to a nat- 
ural frequency of uj n = y/(J 2 /4) — V 2 ~ 0.5 J. Note that 

in this case, owing to the small ratio of Td/J, the varia- 
tion of T f does not affect u n , at least to leading order in 
Tf/J. Thus we expect that the contribution to the DFT 
would arise from the drive frequency u — 0.1J and the 
natural frequency u) n — 0.5J. This explains the double 
peaked structure of the DFT spectrum in panel (a) of 
Fig. Q33 

Finally, let us consider the classical equations at = 
0.55 J. Here the fixed point at (tt/2, 0) which corresponds 
to the energy minima of H (Eq. J5J), is quite shallow 
and even a small transverse drive field is enough to drive 
the system far away from the fixed point. Therefore we 
expect the system to traverse a large part of the phase 
space during its motion even if T a is small. Consequently, 
the DFT spectrum of cos 2 9 is expected to have contribu- 
tions from a wide range of frequencies which translates 
to noisy behavior in the time domain. This expectation 
is verified by numerical computation as can be seen from 
panel (b) of Fig. Note that this noisy nature of the 
AC dynamics of the EOC occurs due to the presence of 
a shallow energy minima; this in turn is the consequence 
of the presence of a QCP at Td/J — 0.5 where the fixed 
point structure of Eq. i|25|) changes. Therefore we expect 
that such a feature will be generic for QCPs. 



FIG. 10: Plot of the DFT power spectrum of the classical 
EOC cos 2 9 for different values of Tj,, obtained by solving the 
classical equations of motion in Eq. (iLTil . Panels (a) and 
(c) show the DFTs in the ferromagnetic and paramagnetic 
phases respectively; panel (b) shows the same for the critical 
region. Note that the involvement of many frequencies in the 
dynamics near the critical point survives even in the classical 
limit. 

Let us first consider the paramagnetic phase where 
Td > J/2 and the fixed point is given by (9f,<j>f) = 
(tt/2, 0). Linearizing Eqs. (|25[> about this fixed point, 
one gets a natural frequency of ui n — yT/ (r^ — J/2) 
for small amplitude motion. Note that this analysis is 
only approximate since for our case Tf itself varies in 
time. But for a small drive amplitude and frequency, 
where T d - T a = 0.8 J < T f < T d + T a = 1.2 J, we expect 
the main contribution to the periodic motion of 9 to come 
from approximately the range 0.49 J < uo n < 0.91J. Thus 
the power spectrum for the motion of (Sf ot ) 2 (t)/S 2 — 
cos 2 9, which has half the time period as cos 9 in the para- 
magnetic phase (since cos 9 runs over both positive and 
negative values), is expected to have significant contri- 
bution from frequencies in the range 0.98 J < u n < 1.8 J 



V. DISCUSSION 

The results obtained in this work can be relevant for 
certain experimental systems. One possible class of sys- 
tems where these results may be applicable are ones with 
long range dipole-dipole interactions such as KH2PO4 or 
Dy(C2HsS04)39H20 2 which exhibit order-disorder tran- 
sitions driven by tunneling fields. However, the effect of 
finite-range of the interactions on the dynamics needs to 
be studied carefully before quantitative predictions can 
be made. We have left this as a subject of future study. 

The other class of systems where our results might 
be applicable are two-component BECs where the inter- 
species interaction (U) is strong compared to the intra- 
species interaction (£/'). Such BECs, in the presence of 
an external RF drive and a fixed chemical potential, can 
be described by an effective Hamiltoniani& 



H, 



err 



x sl + nns x 



(33) 



where the spin operator S z — (&|fo| — b\bi)/2 and S x 



(blb\ + b\bl)/2 are the z and x components of the ef- 
fective spin in terms of the boson operators b^i, | and j 
determines the species or pseudospin index, x = {U — V) 
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is the coefficient of the S"l term, and is the frequency 
of the external RF field. Note that V can be made arbi- 
trarily strong compared to U by tuning the system near 
a Feshbach resonance as discussed for 41 K — 87 Kb system 
in Ref. 0. Alternatively, such a Hamiltonian (Eq. I33|) 
with x < can also be realized following the method dis- 
cussed in Refs. Il6ll7l The quench dynamics described 
in our work can then be realized by a sudden change 
of the frequency of the RF pulse. A more complicated 
RF pulse shape can also lead to an effective sinusoidal 
time-dependent O whereby the AC dynamics can also be 
realized. We also note here there is a vast body of lit- 
erature on the model described by Eq. I33££*^. Most 
of these works study the entanglement property of the 
model for x > 0. In this regime, the model (antiferro- 
magnetic infinite-range Ising model in a transverse field) 
does not have a second-order quantum phase transition 
and has completely different properties than the ferro- 
magnetic model that we study here. 

In conclusion, we have studied the quench and AC dy- 



namics of the EOC of an infinite range ferromagnetic 
Ising model, both classically and quantum mechanically. 
We have shown that both the quench and the AC dynam- 
ics reflect the presence of the quantum critical point of 
the model. In particular, the AC dynamics of the EOC, 
for a small amplitude and moderate frequency transverse 
AC field, exhibits a noisy behavior near the quantum 
critical point which is qualitatively distinct from its pe- 
riodic behavior away from the critical point in either the 
paramagnetic or the ferromagnetic phase. As we have 
shown, this noisy behavior of the AC dynamics of the 
EOC near the QCP follows from quite general arguments, 
and therefore we expect this behavior to qualitatively 
hold for more realistic models exhibiting QCP. Finally, 
we would like to point out that the results presented here 
may be relevant for Ising systems with long range dipolar 
interactions and two-component BECs. 
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for useful correspondence. 
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